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XLII 


ADDITIONAL APPLICATIONS OF THE THEORY OF 
ALGEBRAIC QUATERNIONS TO DYNAMICAL QUESTIONS 


Communicated 8 December 1845. 
[Proc. Roy. Irish Acad. vol. 111 (1847), Appendix, pp. li-lx.] 


The following is the communication made by Sir William R. Hamilton on some additional 
applications of his theory of algebraic quaternions. 

It had been shown to the Academy, at one of their meetings* in the last summer, that the 
differential equations of motion of a system of bodies attracting each other according to 
Newton’s law, might be expressed by the formula: 


dæ m+Am (1) 
dt? ~ —Aa,/(— Aa?) 


in which a is the vector of any one such body, or of any elementary portion of a body, regarded 
as a material point, and referred to an arbitrary origin; m the constant called its mass; « + Aa, 
and m+ Am, the vector and the mass of another point or body of the system; £ the sign of 
summation, relatively to all such other bodies, or attracting elements of the system; and d the 
characteristic of differentiation, performed with respéct to ¢ the time. 

If we confine ourselves for a moment to the consideration of two bodies, m amd m’, and 
suppose r to be the positive number denoting the variable distance between them, so that r is 
the length of the vector æ’ — a, and, therefore, by the principles of this calculus, 


r= a{{—(a' — a)}; (2) 
we shall have, by the formula (1), the two equations, 


da m'r! da’ mr- 


l daa d® wa’ 
which may also be thus written, 
da _ , DE + rda 1 —3 / 
mag mmr (a’—a), m qg mr (a—«a’), 
and which give da dg d?’ dg’ mm’ 
o=m (82 Ga + Gao) +m Q ae tape & tees 


da, da’ being any arbitrary infinitesimal variations of the vectors a, «’, and ôr being the 
corresponding variation of r; because 


cala’ — a) + (x — x) da + dae’ (a — x’) + (— a’) Bx’ 
= — (da — da) (a’ — a) — (a’ — æ) (dae’ — da) = — ô. (x' — a)? =8.7? = 2r dr = — 28. r. 


* [See XLI, equation (A).] 
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And by extending this reasoning to any system of bodies, we deduce from the equation (1) 
this other formula, by which it may be replaced: 


qd? a d?a 
di + de 

Although it is believed that this result (3), if regarded merely as a symbolic form, is new, as 
well as the method by which it has been here obtained; yet if we transform it by the introduc- 
tion of rectangular coordinates, x, y, z, making for this purpose 


4.™ (o da) + “ee ==), (3) 


a=ix+jy+kz, a =ia'+jy' +k’, (4) 


and eliminating the squares and products of the three imaginary units, i, j, k, by the nine 
fundamental relations which were communicated to the Academy in 1843, namely, 


fang? Lae — 1; y=k, jk=i, b=); j= =k, kj=-i, tkh=-3; (5) 
we are conducted, from the equation (3), to a well-known equation of Lagrange, which may 
be written thus: | di d?y i 

Bm (Spa d+ G4 Oy + a); 


where r, by (2), (4), (5), is equal to the known expression, 


r= {{(x' — 2x)? + (y’—y)? + (z —2)9}. 

Ifthe law of the attraction were supposed different from that of the inverse square, a different 
function of r, instead of r-t, should be multiplied by the product of two masses. 

But further, it is not difficult so to operate on the formula (3), as to deduce from it another 
equation which shall be equivalent to the forms that were proposed by the present author, in 
his papers ‘On a General Method in Dynamics’ (published in the Philosophical Transactions) ,* 
as being, at least theoretically, forms for the integrals of the differential equations of motion of 
any system of attracting bodies. For if we observe, that by the principles of the calculus of 
variations, combined with those of the method of vectors, we have the identity, 


da d2a + d2a da = d(da da + da da) — d(da?); 


and if we write v= f -— (=) | 4 (6) 


denoting by v the magnitude or degree (but not the direction) of the velocity of the body of 
which the vector is «; we may transform the fen i into the following: 


d Dhat da 


ymm’ 


which, when operated upon by the characteristic fa dt, that is, when integrated once with 
0 


respect to the time from 0 to t, becomes 


da dy da doy 
x. 3 (82g: ai — 6%)—— apt de da — E da) +dF=0, (8) 

if we make for abridgment „mm 
F= -f di (2z mm) (9) 


* 1834, Part II. 1835, Part I. [See Mathematical Papers, vol. 11, papers II, III.] 
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and denote by da, and a the values which the variation of æ, and the differential coefficient 


of that vector taken with respect to t, are supposed to have at the origin of the time. The 
definite integral denoted here by the letter F is the same which was denoted by the letter S in 
the Essays already referred to, and which was called, in one of those Essays, the Principal 
Function of the motion of a system of bodies; and if we now regard it as a function of the 
time ¢, and of all the final and initial vectors, a, a’ ...,a@ ,a,... of the various bodies of the 
system, and suppose (as we may) that its variation, taken with respect to all those vectors, is 
determined by an equation of the form, 


0=26F + X(cda— oda + da.0 — ôd. To), (10) 


in which o, o, are vectors, we are conducted, by comparison of the coefficients of the arbitrary 
variations of vectors, in the equations (8) and (10), to the two following systems of formulae: 


da > age $ 
rT di =O 5 3 (11) 
day _ rda _ ; j 

ma =% m e ia a's (12) 


of which the former may be regarded as intermediate, and the latter as final integrals of the 
differential equations of motion. The determination of the (vector) coefficient ø, from the 
variation of the (scalar) function F, is an operation of the same kind as the known operation of 
taking a partial differential coefficient, and may, in these new calculations, be called by the 
same name; but in order to be fully understood, it requires some new considerations, of which 
the account must be postponed to another occasion. 

Consider a system of three attracting masses, m; m’, m”, with their corresponding vectors, 
a, a’, «”; and make for abridgment «’—a=/, and «”—a=y; we shall have, by (1), for the 
differential equations of motion of these three masses, referred to an arbitrary origin of vectors, 


the following: - da ait m" 
dP BP) I 
d?(æ + d?(a +2) _ m m” (13) 
~~ BYP)" B-N E- 
d’(a+y)__ m x m ; 
a», Yall =F) = Pr-O BF 
which give, for the internal or relative motions of m’ and m” about m, the equations: 
ap m+m' ml (Ee — y)! y a: ) 
a B= 8) = B=" = fe 


(ek es m le Jeno a 
. dB GA. M=7- BS =P) 

If we suppress the terms multiplied by m” in the first of these equations (14), or the terms 
multiplied by m’ in the second of those equations, we get the differential equation of motion 
of a binary system, under a form, from which it was shewn to the Academy last summer,* 
that the laws of Kepler can be deduced. If we take account of the terms thus suppressed, we 
have, at least in theory, the means of obtaining the perturbations. 


* [See XLI.] 
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Let m be the earth, m’ the moon, m” the sun; then £ and y will be the geocentric vectors of 
the moon and sun; and the laws of the disturbed motion of our satellite will be contained in 
the two equations (14), but especially in the first of these equations. By the principles of the 
present calculus we have the developments, 


(y—P)te=y tty fy tty fy Py +... (15) 
\(-¥?) By+yBb , BA By +yB 
d ES — — = E a AA ie as 
sa i{-(y-6)} (i y? a ae Bye onn bi 


if then we reject the terms of the same order as m”/?y-+, that is, terms depending on the inverse 
fourth power of the distance of the sun from the earth, the disturbing part of the expression for 
the second differential coefficient, taken with respect to the time, of the moon’s geocentric 
vector, will reduce itself in this notation to the following: 


Pa dat le eka —1) — im” ( —vy2)- =i 
EP Ji- (y -p3 J=" YI As NIR a7) 

This symbolic result admits of a simple geometrical interpretation. The symbol yty 
denotes a vector in the plane of the two vectors £ and y, which has the same length as £, and 
is inclined at the same angle to y, but at the other side of that line; so that y bisects the angle 
between 2 and y—!fy. If then we conceive a fictitious moon among the stars, so situated that 
either the sun, or a point opposite to the sun, as seen from the earth, bisects the arc of a great 
circle on the celestial sphere, between the fictitious and the actual moon (the bodies being here 
treated as points); and if we decompose the sun’s disturbing force on the moon into two others, 
directed respectively towards the extremities of that celestial arc which is in this manner 
bisected: one component force, resulting from this decomposition, will be constantly ablati- 
tious, tending directly to increase the distance of the moon from the earth, and bearing to 
the attractive force in the moon’s undisturbed relative orbit, a ratio compounded of the direct 
ratio of half the mass of the sun to the sum of the masses of the earth and moon, and of the 
inverse ratio of the cubes of the distances of the sun and moon from the earth; and the other 
component force, directed towards the fictitious moon, will be exactly triple of the ablatitious force 
thus determined; provided that we still neglect all terms depending on the inverse fourth 
power of the sun’s distance, as we have done in deducing the equation (17), of which the 
theorem here enunciated is an interpretation. A similar result, of course, holds good, for every 
satellite disturbed by the central body of a system. The theorem admits of being proved by 
considerations more elementary, but was suggested to the author by the analysis above 
described; which may be extended, by continuing the developments (15), (16), to the case of 
one planet disturbed by another, and to a more accurate theory of a satellite. 

Without entering into any farther account at present of the attempts which he has made to 
apply the processes and notation of his calculus of quaternions, or method of vectors, to 
questions of physical astronomy, the author wished to state that he had found those processes, 
and that notation, adapt themselves with remarkable facility to questions and results 
respecting Poinsot’s Theory of Mechanical Couples.* A single force, of the ordinary kind, is 
naturally represented by a vector, because it is constructed or represented, in mathematical 
reasoning, by a straight line having direction; but also a cowple, of the kind considered by 
Poinsot, is found, in Sir William Hamilton’s analysis, to admit of being regarded as the vector 
part of the product of two vectors, namely, of those which represent respectively one of the two 


* [Éléments de Statique (Paris, 1804).] 
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forces of the couple, and the straight line drawn to any point of its line of direction from any 
point of the line of direction of the other force. Composition of couples corresponds to addition 
of such vector parts; and the laws of equilibrium of several forces, applied to various points of 
a solid body, are thus included in the two equations, 


Lp=0; X(ab—fa)=0; (18) 


the vector of the point of application being «æ, and the vector representing the force applied 
at that point being 2. The condition of the existence of a single resultant is expressed by the 


formula, DP. U(aP— Ba) + Elaf- Ba). I=. (19) 
Instead of the two equations of equilibrium (18), we may employ the single formula 
x.af=—ce, (20) 


c here denoting a scalar (or real) quantity, which is independent of the origin of vectors, and 
seems to have some title to be called the total tension of the system. 

In mentioning finally some applications of his algebraic method to central surfaces of the 
second order, the author could not but feel that he spoke in the presence of persons, of whom 
several were much better acquainted with the general geometrical properties of those surfaces 
than he could pretend to be. But, while deeply conscious that he had much to learn in this 
department from his brethren of the Dublin School, as well as from mathematicians elsewhere, 
he ventured to hope that the novelty and simplicity of the symbolic forms which he was about 
to submit to their notice might induce some of them to regard the future development of the 
principles of his method as a task not unworthy of their cooperation. He finds, then, that if 
a and p denote two arbitrary but constant vectors, and if p be a variable vector, the equation 
of an ellipsoid with three arbitrary, and, in general, unequal axes, referred to the centre as the 
origin of vectors, may be put under the following form 


(ap + pa)? — (Bp — pp)?=1. (21) 

One of its circumscribing cylinders of revolution is denoted by the equation 
— (Bp — pb =1; - (22) 
the plane of the ellipse of contact by gin pd ahs (23) 

and the system of the two tangent planes parallel hereto, by 

(ap+pa)*=1. (24) 
A hyperboloid of one sheet, touching the same cylinder in the same ellipse, is denoted by 
coc aah (ap + pa)? + (Bp —ph)= —1; (25) 
its asymptotic cone by (ap + pa)?+ (fp—pP)?=0; (26) 


and a hyperboloid of two sheets, with the same asymptotic cone (26), and with the tangent 
planes (24), is represented by the formula 


(ap + px)? + (Bp —pP)P=1. (27) 

By changing p to p — y, in which y is a third arbitrary but constant vector, we introduce an 
arbitrary origin of vectors, or an arbitrary position of the centre of the surface as referred to 
such an origin; and the general problem of determining that individual surface of the second 


order (supposéd to have a centre, until the calculation shall show in any particular question 
that it has none), which shall pass through nine given points, may thus be regarded as equivalent 
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to the problem of finding three constant vectors, «, P, y, which shall, for nine given values of the 
variable vector p, satisfy one equation of the form 


{a(o—y) + (p—y) a}? + {Blp—)—(p—-y) BP = £1; (28) 
with suitable selections of the two ambiguous signs, depending on, and in their turn deter- 
mining, the particular nature of the surface. It is not difficult to transform the equation (28), 
or those which it includes, so as to put in evidence some of the chief properties of surfaces of 
the second order, with respect to their circular sections. 

The recent expressions may be abridged, if we agree to employ the letters S and V as 
characteristics of the operations of taking separately the scalar and vector parts of any 
quaternion to which they are prefixed; for then we shall have 


ap +pa=2.ap, pp—ph=2V.fp; (29) 
so that, by making for abridgment 2% =g&', 28=/’, the equation (21) of the ellipsoid (for 
example) will take the shorter form, 

(S.a'p)?— (V .2'p}=1. (30) 
Another modification of the notation, which, from its geometrical character, will often be 
found useful, or at least illustrative, may be obtained by agreeing to denote by the geometrical 
symbol BA the vector J — a, which is the difference of two other vectors « and # drawn to the 
two points A and B, from any common origin; so that BA is the vector to B from A. Denoting 
also by the symbol CBA the quaternion CB x BA, which is the product of the two vectors 
CB and BA; by DOBA the continued product DC x CB x BA, and so on: the foregoing equa- 
tions of central surfaces may be transformed, and a great number of geometrical processes and 
results expressed under concise and not inelegant forms. For example, the symbols 


V.ABC 


ne To ae (31) 
S. ABCD 
apa V.ABC ia 


will denote, in length and in direction, the perpendiculars let fall, respectively, from the 
summit B on the base AC of a triangle, and from the summit D on the base ABC of a tetra- 
hedron; the sextuple area of this tetrahedron ABCD being expressed in the same notation by 
the symbol 8. ABCD. i 

The developments (15) and (16), with a great number of others, may be included in a formula 
which corresponds to Taylor’s theorem, namely, the following: 

dira? 
fla+da)= (1454355 +...) fas (33) 

the only new circumstance being, that in interpreting or transforming the separate terms, for 
example, the term 4d?fæ, of the resulting development of the function f(æ + dæ), if æ and its 
differential da denote vectors, we must in general employ new rules of differentiation, having 
indeed a very close affinity to the known rules, but modified by the non-commutative character 
of the operation of multiplication in this calculus of vectors or of quaternions. It is thus that, 
instead of writing d.a?=2ada, d.a?=2ada, we have been obliged to write 


d.a*=a.da+da.a; (34) 
0.0%=a.da+da.a. (35) 
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